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1. Introduction 


My review lecture has focused on some recent progress in the physics of brane-worlds. 
Both model building and cosmological implications were considered. Since there is 

Ir eview 

already a published review on such issues l|]J, and since space here is scarce, I have 
decided to develop a particular aspect of brane cosmology, which is both important 
and little studied, namely that of brane-bulk energy exchange. I will review recent 
progress and present some new results. There has been some study of cosmological 
aspects of brane-bulk energy exchange in the recent past Ifl-ljjvJ. Here, I will mostly 

Ijf Kt.t.Z 

focus on the approach described in 

In any theory, where gravity is higher-dimensional classically (around flat space), 
a mechanism must be invoked to explain why observable gravity is four-dimensional 
at experimentally verifiable distances. The easiest and most popular mechanism 
is compactihcation, but alternatives have been also considered lately, namely RS 
localisation (4d in the IR) IQ and brane-induced gravity IQ (4d in the UV). In all 
such realizations, the 4d graviton is accompanied by KK modes that propagate in 
the bulk and can thus mediate energy-exchanging processes between brane and bulk 
matter. Such processes are important because on the one hand they provide stringent 
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constraints on brane-world models and on the other as it was argned in IQ they may 
generate interesting cosmological effects. 

In IQ the cosmological brane-bnlk energy exchange processes , due to KK gravi¬ 
tons were studied in the simplest RS cosmological evolution supplemented with a 
four-dimensional (induced) Einstein term on the brane . It was shown that, in agree¬ 
ment with intuition, such processes do not affect the cosmological evolution during 
4d eras (H ~ p ), while they are important in higher-dimensional eras (H ~ p). 
Moreover, in many important cases, the rate of energy loss can be written as a power 
of the driving brane-energy density, something that can be used in more general 
situations. 

To generalize, we consider more bulk fields and more general bulk-brane cou¬ 
plings. Although this can be formulated by the standard action principle and the 
relevant exact equations derived and studied we will take a short-cut. We will an¬ 
alyze the regime where the bulk energy, at the brane position, can be consistently 
neglected from the equations. Moreover, we parameterize appropriately the energy 
exchange as a specific power of the matter density of the brane. 

The effects found, include accelerated evolution with standard w = 0 matter and 
outflow, inflationary fixed points in the presence of inflow, and tracking behavior 
between observable and mirage (dark) energy. 


2. The model 


We shall be interested in the model described by the action 


S = 




( 2 . 1 ) 


001 


where R is the curvature scalar of the five-dimensional metric Qab-, A,B — 0, 1, 2 , 3, 5, 
A is the bulk cosmological constant, and g a ,g, with a, (3 = 0,1,2,3, is the induced 
metric on the 3-brane. We identify (x,z) with (x, — z), where z = X 5 . However, 
following the conventions of IQ we extend the bulk integration over the entire in¬ 
terval (— 00 , 00 ). The quantity V includes the brane tension as well as quantum 
contributions to the four-dimensional cosmological constant. 

We consider an ansatz for the metric of the form 


ds 2 = —n 2 (t, z)dt 2 + a 2 (f, z)^fijdx l dx j + b 2 (t, z)dz 2 , 


( 2 . 2 ) 


metric 


where 7 ^ is a maximally symmetric 3-dimensional metric. We use k to parameterize 
the spatial curvature. 

The non-zero components of the five-dimensional Einstein tensor are 
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Primes indicate derivatives with respect to z, while dots derivatives with respect to 
t. 

The five-dimensional Einstein equations take the usual form 
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einstein 


where Tac denotes the total energy-momentum tensor. 

Assuming a perfect fluid on the brane and, possibly an additional energy-momentum 
Tc\ m , B in the bulk, we write 


T A C = T A C 


v,b 
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where p and p are the energy density and pressure on the brane, respectively. The 
behavior of T A | mi # is in general complicated in the presence of flows, but we will not 
specify it further at this point. 

We wish to solve the Einstein equations at the location of the brane following 

etruy 

[7 We indicate by the subscript o the value of various quantities on the brane. 


Integrating equations ( |2df ), (pbf|) with respect to z around z = 0 gives the known 
jump conditions 
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The other two Einstein equations (|2.h|) , 
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where To 5 ,T 55 are the 05 and 55 components of T^c I m r evaluated on the brane. 
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Substituting (j2.12|), (i2T13j|) in equations (j2.14|), (Rib) one obtains 
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I 11 the model that reduces to the Randall-Sundrum vacuum IMl in the absence 
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of matter, the first term on the right hand side of equation ( [2. 1 V| ) vanishes. A new 
scale K may be dehned through the relations V = —A /K = 12M 3 K. Sometimes, 
we will relax this condition. 

At this point we find it convenient to employ a coordinate frame in which 
b Q = n Q = 1 in the above equations. This can be achieved by using Gauss nor¬ 
mal coordinates with b(t,z ) = 1, and by going to the temporal gauge on the brane 
with n a — 1. The assumptions for the form of the energy-momentum tensor are then 
specific to this frame. Using [3 = M -6 /144 and 7 = VM~ 6 / 144, and omitting the 
subscript o for convenience in the following, we rewrite equations ( ^iGj ) and 
in the equivalent first order form 


p + 3(1 + w) - p = -T 
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= (3p 2 + 27 (p + x) ~ 4 + A 
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where p = wp, T = 2T ° 5 is the discontinuity of the zero-five component of the bulk 
energy-momentum tensor, and A = (A + U 2 / 12 M 3 )/ 12 M 3 the effective cosmological 
constant on the brane. 
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In the equations above, Eq. (pT18|) is the definition of the auxiliary density 


With this definition, the other two equations are equivalent to the original system 

IIP .1 11 ’A ■ I 


7T7|). As we will see later on, in the special case of no-exchange (T = 0) x 


represents the mirage radiation reflecting the non-zero Weyl tensor of the bulk. 
The second order equation ( jklY| ) for the scale factor becomes 


- = -(2 + 3 w)/3p 2 - (1 + 3w)xp - 2yy + A. 
a 


( 2 . 20 ) 


As mentioned above, in the Randall-Sundrum model the effective cosmological con¬ 
stant A vanishes, and this is the value we shall assume for most of the rest. 
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We may now pass to dimensionless quantities. The AdS scale is defined as 
V = M 3 K and the dimensionless densities as 


P 72M 3 K ’ X 72 M 3 K 


A A | 1 

K 2 ~ 12 M 3 K 2 + 144 
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as well as a rescaled time r = K t. The combination k = kjK 2 is then dimensionless 
and by scaling a we can set it to k = 0, ±1. Then the cosmological equations become 
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p + 3(1 + w) -p = -T , — = 36 p~ + p + X -y + A 
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X + 4-x = (72p+ 1)T--T 5 , , q = - = —36(2 + 3w)p 2 --— P~X~^T b + A 
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(2.24) 

where now dots stand for derivatives with respect to r. 

We may separate the dynamics associated with T ° 5 and T 5 5 by introducing a 
new energy density a so that 
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In order to derive a solution that is largely independent of the bulk dynamics, the 

T\ term on the right hand side of the same equation must be negligible relative to the 

second one. This is possible if we assume that the diagonal elements of the various 
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contributions to the energy-momentum tensor satisfy the schematic inequality 
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Our assumption is that the bulk matter relative to the bulk vacuum energy is much 
less important than the brane matter relative to the brane vacuum energy. In this 
case the bulk is largely unperturbed by the exchange of energy with the brane. 
When the off-diagonal term T ° 5 is of the same order of magnitude or smaller than 

tor. Horn 

the diagonal ones, the inequality (|2.26|) implies T pK. 
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Thus, under this assumption, we may neglect T 5 from ( |2724| ). Moreover, we 
will parameterize T as a function of the driving energy density p. The precise form 
depends on the details of the interaction between brane and bulk. In the AdS scaling 
region this has typically a power dependence T ~ p v . 


2.1 Other forms of the cosmological equations 
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Combining equations ( |2.28| ) we obtain 


a ~T = _ 3(1 + w)p - eT(p) (36 p 2 + p + X~ -5 + A 
da \ a z 
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Similarly, equations ( )2 23j ) and ( P^4| ) give 
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where e = 1 refers to expansion, while e = —1 to contraction. These two equations 
form a two-dimensional dynamical system. The function y(p) is obtained from the 
equation 


3(1 + w)p\j 36 p 2 + p + y- -^ + A + eT(p )j ^ 


= 4\J 36p 2 + p + x- 4 + A- e (36p + 1) T(p). 


(2.29) 


Note that the equations of contraction are those of expansion with the roles of outflow 
and influx interchanged. 

For k = 0, T = p~ 3 / 2 T and ( = ,/36p + - + 1 + - the equation can be simplified 
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Finally the acceleration satisfies 
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2.2 The four-dimensional regime 


When p « 1 we are in the 4d regime ( H 2 ~ p ). Here the cosmological equations 
simplify 


For k 
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3. Inflating fixed points 


An interesting feature of the cosmological equations is the possible presence of ac¬ 
celerating cosmological solutions. We may look for exponential expansion with a 
constant Hubble parameter H, even if the brane content is not pure vacuum energy. 

We will restrict ourselves for simplicity to the 4d regime. For the non-linear analysis 
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we refer the reader to 


For a fixed point equations (|2 ~3 


without necessarily requiring w = — 1 (A = k = 0) 
(denoted by *) of these equations for k = 0 satisfy 


75ff) must have a time-independent solution, 
The possible fixed points 


3//.(l +w)p, = -T(p,) , II* = [>, + x, , 4 H'X*=T(p,)- 


(3.1) 



It is clear from equation (^T|) that, for positive matter density on the brane (p > 0), 
flow of energy into the brane (T(p) < 0) is necessary. 

The accretion of energy from the bulk depends on the dynamical mechanism that 
localizes particles on the brane. Its details are outside the scope of our discussion. 
However, it is not difficult to imagine scenaria that would lead to accretion. If the 
brane initially has very low energy density, energy can by transferred onto it by bulk 
particles such as gravitons. An equilibrium is expected to set in if the brane energy 
density reaches a limiting value. As a result, a physically motivated behavior for the 
function T(p) is to be negative for small p and cross zero towards positive values for 
larger densities. In the case of accretion it is also natural to expect that the energy 
transfer approaches a negative constant value for p —> 0. 
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The solution of equations 


satisfies 


T (p*) = ~2 ( X + w )^ 1 ~ 3w P * f2 
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For a general form of T(p) equation ( p+j ) is an algebraic equation with a discrete 
number of roots. For any value of m in the region — 1 < w < 1/3 a solution is possible. 
The corresponding cosmological model has a scale factor that grows exponentially 
with time. The energy density on the brane remains constant due to the energy flow 
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from the bulk. This is very similar to the steady state model of cosmology ITT]. The 
main differences are that the energy density is not spontaneously generated, and the 
Hubble parameter receives an additional contribution from the “mirage” field x (see 
equation (^ 1 |)). 

The stability of the fixed point ( |^1| ) determines whether the exponentially ex¬ 
panding solution is an attractor of neighboring cosmological flows. If we consider a 
small homogeneous perturbation around the fixed point (p = p* + 5p, x — X* + ^x) 
we find that 5p, 5\ satisfy 
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and we have employed the relations ( 
matrix A4 are 


and T(p) oc p v . The eigenvalues of the 
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7 + 3m — 3z/(l + m) ± 3 + 2 z/)(l + m) + [7 + 3m — 3z/(l + w)Y 
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(3.6) 

For —1 < m < 1/3, 0 < v < 3/2 they both have a positive real part. As we 
have assumed T(p) < 0, the fixed point is stable in this case. The approach to the 
fixed-point values depends on the sign of the quantity under the square root. If 
this is negative the energy density oscillates with diminishing amplitude around its 
fixed-point value. 
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4. Tracking solutions 

We will now analyze the case v — 3/2 which lies at the boundary of the stability 
region discussed above. We will thus assume that T = A p 3 / 2 ; and that the universe 




















expands and is dominated by non-relativistic matter (w=0). 


equation 


becomes (A = 0) 


Then, in the 4d regime, 
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We will parameterize the dimension-less coefficient A a,s A = p — A is de¬ 
termined by the details of the microscopic cross section that gives rise to this type 
of energy exchange, p running on non-negative real numbers parameterizes all pos¬ 
sible values of A. When we have expansion, A > 0 means outflow. When we have 
contraction A < 0 means outflow. 

[track 

The general solution of equation (|f.l|) is 
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Since H 2 = p + x the equation above can be re-written, in terms of £ = H/ Wo. 

e V 
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and can be integrated as a function of a with the result 
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p(a ) can be obtained by solving (p^|) and substituting into (pf|) ■ Finally \{a) = 


P(C 2 -1). 


We will first study a few special cases : 
(i) p = 1. Here we obtain 


“ 1 + Cp l/3 => H 2 — p + C pt 3 , , X = ~ (4.5) 

compatible with the absence of energy exchange in this case and consequent inde¬ 
pendence of the evolution of p, x- 
(ii) /i = -1/2. 

C = 2 + Cp^ H 2 = p{2 + Cp 1/6 ) 2 (4.6) 

For large a 

p= {iSf 2 a_l5/4+ ° ( a_l5/2 ) (47) 

There is a similar solution for p = 1/2 with asymptotic contraction. The cases 
p = ±2 are equivalent to the above. 
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(iii) 1/2 < /i. It corresponds to — | < A. Asymptotically (a —> oo) we obtain 
the tracking solution 






( = p + C + ... , H 2 = /Ap + 2pC p2(4 M ^'-i) +1 + ... (4.8) 

p ~ C" a^ -4 + ... , x = {p 2 - l)p + ... (4.9) 

Here, although the initial conditions for the real p and mirage x energy density are 
arbitrarily different (parameterized by the independent integration constants C,C), 
at late times the scale similarly with the scale factor 
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p~C' + , ±=(,i 2 _!) + .. 
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Tims, the dark energy behaves as the visible energy, and such a mechanism could be 
used so that bulk energy simulates dark matter. 

(iv) — 2 < p < 0 It corresponds to — § < A. For asymptotically small p we 


obtain 
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C = ^T + Cp^ z ^ + ... , H' 1 = ^ + 2f- P^^ +1 + ... 

\p\ p 2 \p\ 


On the other hand 


p ~ C 1 a" 1 4 + ... , x = ~ T~ P+ ..• 

fi- 
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Here we have again tracking behavior. 

All other ranges have asymptotic ( which is negative and thus unphysical. 
Finally, in the case of outflow, there is a fixed point in the 5d regime, when 
A 2 < 9/4 with 

9-4A 2 „ A 1 9 - 4 A 2 

P * 8.81 ’ * 108 V 2 ' ' ^ 

This is a saddle point 

It is interesting to note that a similar_tracking behavior has been observed in 
matter interacting with the dilaton in 


vaQ 

ez 

12 



, (§%]]) with H = FT* 


5. Fixed points in the non-linear regime 

We will consider solutions to the non-linear system 
constant. The equations also imply that p = p*, T — T*, x — X* are a l so constant. 
We will see that although there may be a leftover cosmological constant A on the 
brane, the cosmological acceleration because of energy inflow, may be much smaller 
than \f\. 

From the equations we obtain 
1 


pf = 


144(1 + 3w ) 


(1 - 3w ) ± y/(l - 3w) 2 + 1152(1 + 3m)(A - H 2 ) 


(5.1) 


10 






















X: 


* = —^(1 + W )P* [? 2 p* + 1] , T* = -3(1 + w)H*p* , q* = Hi (5.2) 

Assuming the rate of expansion H* to be small compared to the cosmological 
constant A, we have the following two possibilities 

(i) A dominates in the square root. In this case p* ~ ^ i8(i+3n>j ~ There is still 
space for this approximation to be correct and p* << 1 so that we are in the 4d 
period. 

(ii) In the opposite case p* ~ 72(1+3^) an d we can t> e either i n the 4d or the 5d 
regime. 

I11 either case , energy exchange can mask a leftover brane cosmological constant 

6. Other accelerating solutions 

We will present here two different families of solutions that are characteristic in their 
classes. 


a 

a 



Figure 1: Outflow, k = 0, w = 0, v = 1. The arrows show the direction of increasing 
scale factor 


A global phase portrait of q = a/a with respect to p during expansion in the 
outflow case for k = 0, w = 0, v = 1 is shown in Figure 1. All solutions are below 
the limiting parabola q < ^P - 36(1 + 3 w)p 2 . 

One recognizes two families of solutions: The first have q < 0 for all values of p, 
while the second start with a deceleration era for large p, enter an acceleration era 
and then return to deceleration for small enough values of p. 

Solutions corresponding to initial conditions with positive q (always under the 
limiting parabola shown with the dotted line), necessarily had a deceleration era in 
the past, and are going to end with an eternal deceleration era also. The straight 
dashed line represents the standard FRW solution without the effects of energy ex¬ 
change. 
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The global phase portrait of g = a/a with respect to p during expansion for 
the case k = 0, w = 0, v = 1 is shown in Figure 2. The presence of the limiting 
parabola as in the outflow case is apparent. However, new characteristics appear. 
For example, pi ^ attracts to eternal acceleration a whole family of solutions which 
start their evolution at either very low or very high densities. There is another family 
of solutions which are attracted to acceleration by p * and which eventually exit to 
a deceleration era. Finally, there is a family of solutions, near the limiting parabola, 
which start with acceleration at very low densities, and eventually exit to eternal 
deceleration, while their density increases monotonically with time because of the 
influx. 



For v 1 one expects a different set of fixed points with varying behaviors 
around them. 

7. Conclusions 

An approximate phenomenological analysis of the cosmological effects of brane-bulk 
energy exchange indicates the presence of several interesting phenomena ranging 
from acceleration (and its exit) from energy outflow, to inflationary fixed points in 
the case of inflow to tracking solutions. 

This makes imperative the detailed study of the associated cosmological dy¬ 
namics. I particular the role of the ’’small-back-reaction approximation” must be 
understood. It is expected that this is ok at late times i the evolution but it breaking 
at large densities should be quantified. 

A particularly interesting direction is the holographic formulation of this prob¬ 
lem. This should correspond to the observable matter, interacting directly to classical 
four-dimensional gravity as well as to a ’’hidden” four-dimensional gauge theory (a 
perturbation of N=4 super Yang-Mills). Brane-bulk energy exchange corresponds 
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to the interaction between observable matter and ” hidden” matter. A quantitative 
study should clarify several aspects of the problem. 
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